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1  Problem  Studied 


The  structure  of  the  eddies  in  the  atmospheric  boundary  layer  (ABL)  that  develops  above  heterogeneous 
land  surfaces  seems  to  be  significantly  affected  by  the  spatial  distribution  of  the  land-surface  sensible 
heat  flux,  and  by  topography.  This  is  expected  to  have  an  important  impact  on  the  dispersion  of 
materials  (i.e.,  gases,  aerosols,  and  other  particles)  in  the  ABL. 

The  research  project  conducted  liere  consisted  of  three  tasks,  which  were  performed  with  the  Col¬ 
orado  State  University  -  Regional  Atmospheric  Modeling  System  (RAMS),  a  state-of-the-art  modeling 
system  designed  to  support  basic  research  in  the  atmospheric  sciences:  (i)  An  evaluation  of  the  Large- 
Eddy  Simulation  (LES)  option  of  RA.MS  (hereafter  referred  as  "RAMS-LES”)  using  the  the  First  In¬ 
ternational  Satellite  Land  Surfat  e  Climatology'  Project  (ISLSCP)  Field  Experiment  (FIFE)  data  sets; 
(ii)  A  systematic  analysis  of  the  impac  t  of  microscale  and  meso-7-scale  heat  patches  and  topographical 
features  on  the  structure  of  the  API,  w  ith  the  RAMS-LES;  and  (iii)  An  analysis  of  the  impact  of  these 
landscape  heterogeneities  on  the  di.sperston  of  materials  in  the  ABL.  This  third  task  was  performed 
with  the  Lagrangian  Particle  Dispersum  .Model  (LPDM)  option  of  RAMS  “fed”  with  the  atmospheric 
conditions  generated  by  the  large-eddv  .simulations  performed  under  Task  1  and  Task  2. 

This  study  contributes  to  a  bet  tin  understanding  of  the  physical  mechanisms  involved  in  the  disper¬ 
sion  of  materials  in  the  ABL  that  (lo\elops  above  heterogeneous  (realistic)  land  surfaces.  In  addition, 
the  impact  of  land-surface  heterogeneity  on  the  mean  atmospheric  variables  in  the  ABL  provides  new 
insights  on  the  representation  of  the  .ABL  in  atmospheric  models  which  do  not  have  the  very-high  grid 
resolution  required  to  resolve  turbulent  eddies. 

2  Summary  of  the  Most  Important  Results 

2.1  Model  Validation 

-As  part  of  Task  1,  the  RAMS-LES  was  used  to  simulate  the  convective  boundary  layer  (CBL)  that 
developed  on  July  1.  1987,  over  the  FIFE  domain.  Three  simulations  were  produced  using  different 
boundary  conditions  at  the  ground  surface,  namely;  (i)  spatial  distribution  of  topography  and  spatial 
distribution  of  surface  heat  fluxes:  (ii)  spatial  distribution  of  topography  but  mean  surface  heat  fluxes; 
and  (iii)  no  topography  and  mean  surface  heat  fluxes.  The  diurnal  variation  of  mean  surface  fluxes  and 
their  spatial  distribution  were  derived  from  the  FIFE  network  of  observations.  In  all  cases,  the  model 
was  initialized  with  the  atmospheric  sounding  observed  in  this  domain  at  7  a.m.,  and  run  until  3  p.m. 
The  resulting  mean  profiles  of  temperature  and  specific  humidity  were  compared  to  those  observed  with 
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atmospheric  soundings  at  9  a.m..  10:30  a.m.,  and  12:30  p.m.  The  simulated  structure  of  turbulence 
was  qualitatively  compared  with  that  obtained  from  a  volume-imaging  lidar  (VIL)  scanning  the  CBL 
over  the  simulated  domain,  during  that  day.  Power  spectra  and  autocorrelations  of  mixing  ratio  were 
calculated  from  the  model  outputs  and  were  compared  to  those  obtained  from  the  VIL. 

Overall,  the  model  performed  quite  well.  Observed  atmospheric  soundings  were  within  1  K  and 
1  g/kg  of  the  simulated  mean  profiles  of  temperature  and  specific  humidity,  respectively,  and  indicated 
that  the  model  correctly  predicts  the  CBL  height.  Similarities  in  the  structure  of  the  eddies  obtained 
from  the  model  and  the  VIL  were  dearly  identified.  Spectral  analysis  indicated  that  resolved  eddies 
(he.,  eddies  larger  than  200  m)  are  relatively  well  simulated  with  the  model,  but  that  the  energy 
cascade  is  not  well  represented  1)\’  the  Deardorff  1.5-order-of-closure  subgrid-scale  parameterization. 
Autocorrelation  analysis  indicated  that  the  model  correctly  simulates  the  characteristic  size  of  the 
eddies,  but  that  their  mean  iifetnne  is  longer  than  that  observed  with  the  VIL,  indicating  a  too  weak 
dissipation  of  the  eddies  by  the  sui)2,rid-scale  scheme.  Thus,  this  study  emphasized  the  needs  to  develop 
better  subgrid-scale  parameterizatiuns  for  LES  models.  The  different  simulations  also  indicated  that 
topographical  features  of  tlie  cjrder  of  100  m  and  micro-d  scale  heterogeneity  of  surface  heat  fluxes  had 
only  a  minor  to  modest  impact  on  the  CBL  developing  over  a  relatively  humid  surface. 

A  complete  description  of  this  evaluation  study,  including  a  discussion  of  the  dataset  produced  by 
the  VIL.  is  given  in  Avissar  rt  nl.i  1998).  This  paper  is  provided  in  Appendix  A. 

2.2  Impact  of  Landscape  Heterogeneity 

In  addition  to  the  evaluation  of  RAMS-LES,  the  above-mentioned  numerical  experiments  indicated  that 
topographical  features  of  the  order  of  100  m  and  micro-d  scale  heterogeneity  of  surface  heat  fluxes  had 
only  a  minor  impact  on  the  CBL  developing  over  a  relatively  humid  surface  (see  Avissar  et  a/.,  1998). 
While  these  simulations  represented  real  topographical  features,  they  provide  only  limited  insights  on 
the  key  issue  being  investigated  here,  namely  the  impact  of  landscape  heterogeneity  on  the  ABL. 

As  part  of  Task  2,  the  effects  on  the  CBL  of  surface  heterogeneities  produced  by  various  topograph¬ 
ical  features  as  well  as  sensible  heat  waves  with  different  means,  amplitudes,  and  wavelengths,  were 
investigated.  We  found  that  the  impact  of  topography  and  of  heat-wave  amplitude  and  wavelength  is 
nonlinearly  dependent  upon  the  mean  heating  rate.  The  circulations  (or  rolls)  resulting  from  surface 
heterogeneity  are  relatively  strong  when  the  hills  dimensions  (z.e.,  the  characteristic  height  and  width 
of  the  hills)  and  the  heat-wave  amplitude  and  wavelength  are  relatively  large,  especially  at  low  mean 
heating  rate.  In  that  case  the  profiles  of  horizontally-averaged  variables  are  quite  strongly  modified 
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in  the  CBL.  The  potential  temperature  is  not  constant  with  elevation,  and  the  sensible  heat  flux  con¬ 
siderably  departs  from  the  linear  variation  with  height  obtained  in  a  typical  CBL  that  develops  over 
an  homogeneous  domain.  The  mean  turbulence  kinetic  energy  (TKE)  profile  depicts  two  maxima,  one 
near  the  ground  surface  and  one  near  the  top  of  the  CBL,  corresponding  to  the  strong  horizontal  flow 
that  develops  near  the  ground  surface  and  the  return  flow  at  the  top  of  the  CBL. 

This  study  is  described  in  detail  in  Avissar  and  Schmidt  (1998)  and  Gopalakrishnan  et  aL  (2000a), 
which  are  provided  in  Appendix  B  and  Appendix  C,  respectfully.  A  third  paper,  by  Baidya  Roy  and 
Avissar  (2000)  resulting  from  this  task  summarizes  a  wavelet  analysis  of  the  atmospheric  scales  resulting 
from  the  various  characteristic  scales  of  landscape  heterogeneity  forcing.  This  paper  is  provided  in 
Appendix  D. 

2.3  Impact  on  Dispersion  of  Materials 

The  velocity  components  resulting;  from  the  various  runs  produced  as  part  of  Task  1  and  Task  2  were 
used  to  “feed”  the  Lagrangiaii  Particle  Dispersion  Model  (LPDM)  option  of  RAMS,  to  study  the 
effects  of  landscape  heterogeiieit\'  on  the  dispersion  of  materials  in  the  CBL.  The  major  findings  from 
the  various  experiments  performed  in  this  study  are: 

1.  The  turbulent  eddies  that  develop  above  a  flat  ground  surface  heated  homogeneously  are  randomly 
distributed  within  the  CBL.  Neutrally  buoyant  materials  emitted  at  ground  level  of  such  a  domain 
remain  there  momentarily,  hut  then  quickly  “lift-off”  towards  the  middle  of  the  CBL.  When 
released  from  elevated  sources,  they  first  descend  towards  the  ground  surface  and  then  move 
towards  the  middle  of  the  CBL.  While  near-surface  plumes  ascend  due  to  the  “sweep  out”  of 
particles  converging  into  strong  updrafts  generated  near  the  surface,  elevated  plumes  are  most 
likely  affected  by  downdrafts,  whose  area  coverage  is  much  larger  than  updrafts.  These  results  are 
not  new  (see  e.g,,  Lamb.  1984:  Deardorff  1987;  Briggs,  1993),  and  they  are  used  here  to  confirm 
the  correct  behavior  of  the  model  and,  mostly,  as  a  reference  for  the  dispersion  simulations  above 
heterogeneous  domains: 

2.  Heat  flux  heterogeneities  with  a  characteristic  length  scale  larger  than  about  5  km  have  a  remark¬ 
able  influence  on  particle  dispersion  in  the  CBL.  The  horizontal  pressure  gradients  created  by 
these  heterogeneities  impede  vertical  mixing.  For  a  near-surface  release,  particles  are  advected 
horizontally  rather  than  **lifting-off.”  increasing  the  concentration  there.  Particles  released  at 
higher  elevations  reach  the  surface  much  more  slowly  than  when  released  above  a  flat,  homoge- 
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neous  domain.  The  weaker  the  mean  sensible  heat  flux  injected  in  the  atmosphere  at  the  ground 
surface,  the  more  significant  is  this  effect: 

3.  Hilly  terrain  has  only  a  modest  impact  on  the  vertical  mixing  of  passive  materials  in  the  CBL. 
This  is  probably  because  horizontal  pressure  gradients  resulting  from  small  topographical  features 
are  weak  and,  consequently,  buoyancy  appears  to  have  a  stronger  control  on  the  dispersion  of 
particles  released  at  lower  elevations.  This  is  true  even  with  hills  as  high  as  25%  the  height  of 
the  CBL.  However,  big  hills  have  a  noticeable  effect  on  the  dispersion  of  particles  released  from 
higher  elevations.  In  particular,  for  a  source  height  located  at  about  25%  the  height  of  the  CBL, 
the  locus  of  the  maximum  concentration  in  crosswind-integrated  plume  descends  to  the  surface 
of  the  hill  noticeably  more  slowly  than  above  a  flat,  homogeneous  domain. 

This  study  is  described  in  detail  in  Gopalakrishnan  et  al.  (2000b),  which  is  provided  in  Appendix  E. 
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